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Abstract 

In the present work, we study qualitatively the physics of granular temperature 
coexistence, for a mixture of two different species. Our model captures its essential 
aspects and this allows us to get insights on the physical mechanisms of distinct 
temperature coexistence, in a way which is not obscured by the complexities of 
kinetic theories or numerical simulations. Our simple model is consistent with limit 
situations where we should expect equality for the granular temperatures for the 
mixture. 
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1 Introduction 

Granular Systems (GS) are ubiquitous in Nature [1]. They consist of large 
numbers of distinct, inelastic, rough grains. When studying GS, we are faced 
with mesoscopic collective effects, such as frictional forces, that are not present 
in molecular gases [1-3]. The presence of energy dissipation during internal 
collisions [4] make GS very different from molecular gases: true equilibrium is 
only possible when the total kinetic energy is completely dissipated. Thus, in 
order to keep a GS in a steady-state, we need to inject energy through the 
boundaries or by means of a coupled heath bath [5-8]. The parameter used 
to account for the kinetic energy present in a GS is the so called granular 

1 Partially supported by FAPERJ and CNPq, Brazil. 



Preprint submitted to Elsevier Preprint 



1 February 2008 



temperature, i.e., the average kinetic energy per grain. It is clearly related to 
the usual microscopic definition of temperature in statistical mechanics [9]. 

Thus, we define macroscopic granular quantities taking inspiration from ther- 
modynamic concepts familiar to us. For instance, pressure and granular num- 
ber density, among others, are widely used when studying GS [9,10]. However, 
the naive application of thermodynamic concepts to GS can sometimes be 
misleading. Take pressure for instance: for molecular gases, as we increase the 
number density, at fixed energy feeding rate, the pressure should increase. 
However, for GS, the reverse happens for high enough densities [7]. This tells 
us to be careful when applying thermodynamic concepts, such as the granular 
temperature T g , to GS. 

Thermodynamics is based on a few simple principles, much in the same way 
as Euclidean geometry is based on five postulates. However, we can build non- 
Euclidean geometries from the first four postulates and the negation of the 
fifth. We ask ourselves: is there a similar hope for GS thermodynamics? In this 
spirit, we start by studying the most basic law of granular thermodynamics: 
the zeroth law, or the temperature equilibrium law. Our model is a very simple 
one, based on a Langevin-like approach to the collisions of different granular 
species in a mixture. 

We study the transient behavior of a binary granular gas mixture evolving 
toward a steady-state [11]. The mixture components are smooth, inelastic 
identical spheres of two type: a species of mass M and another species of 
mass m < M and same radius. A convenient Langevin-like collision term is 
introduced, aiming to take into account the more complicated aspects of the 
physics in a simplified, effective way. Experimental realizations of granular 
temperature coexistence [12] have shown some characteristics that can be 
reproduced by the present model. In our opinion, the model's main advantage 
is that it permits us to visualize some important physical aspects of an actual 
system in a mathematical setting which is much simpler then that of more 
complete theories [13]. 

This paper is organized as follows. In Section 2, we study the behavior of 
an intruder in a bath of particles. In Section 3, we analyze the mixture of 
granular species while we calculate the temperature ratio in Section 4. We 
conclude briefly in Section 5. 



2 Test particle in a granular bath 

We study the behavior of a ID test particle of mass M, velocity V, immersed 
in granular bath of particles of mass m = 7M and same radius of the test 
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particle [14] . The granular bath interacts with the test particle through frontal 
inelastic collisions characterized by a coefficient of restitution < e < 1. 

We assume that the typical bath particle has a constant driving velocity 

vd and we suppose that collisions happen with equal probabilities for the 
direction of the bath particle's velocity. The choice above is a reasonable one, 
since there are always four possible ways for a bath particle to collide with the 
test particle (with velocity V): two frontal collisions (with V and ±Vd) and 
two rear ones (again with V and ±t>£>). For each case, regardless of whether 
| V| > v D or | V| < v D , a simple analysis shows that the next collision will occur 
either with a +vd or with a —vd bath particle. A further implicit simplifying 
hypothesis is that we ignore the dependence of the collision rate on the relative 
velocity. This is similar to the so called Maxwell model [15]. We thus write 
the momentum conservation and the inelastic behavior equation as 



MV ±mv D = MV' + mv', (1) 
V'-v' = -e(V-(±v D )). (2) 

After a little manipulation, we obtain 



The test particle's change in velocity is thus 

AV = V'-V = /3[-V + rjv D ], where f3 = 7 ^ + ^ , (4) 



and 7] is a random variable with values ±1 with equal probability distribution, 
< >= 0, ?7 2 = 1. 

The test particle's collisional history is determined by its sequence of collisions 
with the bath particles. We assign an independent random variable ?7jv-i for 
the iV-th collision. Then, we can write the velocity of the test particle, after 
the N— th collision as 



JV-l 

V N = (1- (3fV + (3v D £(!- Pf- 1 ^, (5) 

i=0 

where Vo is the initial velocity of the test particle and 

< rji >= 0, < rjiTjj >= Sij, rjl = I. (6) 
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It can be seen from above that the memory of the initial velocity is lost 
exponentially fast. 

Next, we calculate the cases of the test particle being a bath particle (7=1) 
and that of it being a massive intruder (7 7^ 1) [14,16]. In these two limiting 
cases, the quadratic velocities will be in equilibrium with the bath and their ra- 
tio should reflect the ratio between test particle-bath particle quadratic ratio. 
Obviously, we have Vbath 7^ vd since we allow Vb at h to fluctuate while supposing 
vd fixed. Thus, it is straightforward to show that 



vLh = „ lim <Vn> 

iv -^00,7=1 



J test 



lim <V%>: 

N—>co 



D 



7(1 +g) 

2 + 7 (l-e) 



(7) 
(8) 



From Eq. 7-8, we obtain the granular temperature ratio in the steady state 
for a test particle immersed in a granular bath: 

= \m < v 2 bath > = = 2 + 7(1-^) (Q) 

Tf* \M<vl st > fvl^±^ 3-e ' {) 



In order to compare the result of Eq. 9 with the literature on the tracer limit, 
we see that in the quasi-elastic case, e — > 1, we have up to first order in the 
limit of the heavy tracer 7-^0 

rptest ^ 

( 10 ) 



From reference [13], which agrees with references [14,16], we observe that from 
Eq. 31 and Fig. 3, the same limit above is obtained. The present model seems 
to have better agreement with kinetic theories on this region. The case of larger 
7 does not agree so well. However, its aim is not to obtain accurate quantita- 
tive results, but to get insights into the mechanism of granular temperature 
coexistence. 

The cases of: same species (7 = 1) and elastic motion (e — 1) imply that 
T bath = T test^ ag expected. 

We observe that the above difference in granular temperatures is due to the 
factor (3 being smaller for the massive intruder, 7 < 1 [14] , leading to a smaller 
loss of speed during an inelastic collision giving 

(i-A-) = (i-I^^)>(i-A-) = (i-^). 
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(see Eq. 5). This reflects the fact that, due to its inertia, a massive inelastic 
intruder will lose comparatively less energy to a lighter bath than a less massive 
particle would. 

Our next step will be to characterize the mixture of different species. 



3 Granular mixture 



We suppose that the granular fluid under study is composed of homogeneous 
particles of identical diameters. Similarly to the previous section, any chosen 
particle's collisions will occur with a bath of particles. But now this bath con- 
sists of a mix of two types of particles with different constant driving velocities, 
namely vp for particles m and for particles M [11] (the above mentioned 
particle belongs to one of these families). The coefficients of restitution are Ei, 
e 2 and e' corresponding to the mm, MM, mM (or Mm) collision possibilities. 
In what follows, we use the definitions 

1+fi l + £ 2 R 7(1 + , nu 

«i = — , «2 = —3—. P = 1+ = 7C- (11) 



We need two families of random variables in order to describe the chosen 
particle's N— th collision of a: a variable 9 m ^ or 9m,n, to describe with what 
species the collision occurs; another, rf m> N or t]m,n, to give the direction of the 
bath particle. These variables have the following distributions: 9 m ^ = 1 with 
probability x, 9 mt N = with probability 1 — x, 9m,n = 1 with probability 
1 — x, 9m,n = with probability x, where 

x = ^^, (12) 
Pm + Pm 



and p m and pu are the species number densities; the r\— variables are dis- 
tributed with equal probabilities. 

According to the distributions above and skipping the algebra, in the limit 
N — > 00 we obtain: < Voo >= and 



,2 ^ _ xol\vI + (1 - z)C 2 Vg n - 

00 l-(x(l- ai ) 2 + (l-*)(l-C) 2 )' [6> 



Similarly, we have < V^, >= and 
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(I - x)a 2 V 2 + x(3 2 v 2 
> ^ - 1 _ ((l- x )(l-a 2 ) 2 + x(l - /3) 2 ^ 



<c>= , 77^ r/r ^ ^ - (14) 



4 Temperature ratio: general situation 



As in the case of a test particle, we define the granular temperature ratio as 

^ \M<V*> 1 <V^> 1 V 2 D [b) 

where we assumed that the same relationship holds both for the steady-state 
velocities and the driving velocities. 

Replacing Eq. 13 and Eq. 14 into Eq. 15, and solving for yields 

s 



T™ f [xa 2 - (1 - x)a 2 2 A] 



T g M 7 [ 2x/3 2 A 

^[xa\ - (1 - a;)a|A] 2 + 4x(l - x)(3 2 C 2 A 

2x[3 2 A j ' 

where 

A _ l-x(l- ai ) 2 -(l-x)(l-0 2 
1- (l-x)(l-a 2 ) 2 -x(l-/?) 2 



+ 



(16) 



(17) 



We notice, in Eq. 16, that the dependence of 7fw on 7 is stronger than on 

9 

the inelasticity parameters, as mentioned in reference [12]. Some important 
physical limits are listed below. 

Elastic limit. In this limit, e± = e 2 = s' = 1, and we have 

2t T rn 

ai = a 2 = 1, = 7 C = T-f - =► =fc = 1, (18) 

1 + 7 / M 



for any value of 7. The elastic equilibrium regime is such that all granular 
temperatures are the same, as expected. 

Same species, same inelasticity. In this case we set 7 = 1 and e\ = e 2 = s' = £■ 
We have now 

ai = ot 2 = p = C = — 7: — • 
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It is also trivial to check that ^ = 1- But, for 7 7^ 1 and e 1 ,e 2) £' < 1 the 

a 

granular temperature ratio is non trivial, assuming values different from 1. 

Tracer limit. It is straightforward to show that we recover Eq. 9 when x — > 
[14,16]: 

lim g^ = 2 + 7 (l- £ ) 



Comparison with experiments. In a recent work, Feitosa and Menon [12] used 
beads of equal diameters, made from different materials in order to study the 
granular temperature behavior of a granular mix. According to their experi- 
mental data, for a 1/2-1/2 mix of glass-brass: x = 0.5, 7 = 0.28, £ 1 = 0.83, 
e 2 = 0.61, e' = 0.72, ai = 0.92, a 2 = 0.81, f3 = 0.38 and C = 1-34. We assumed 
e' = £1 + £2 in order to carry out the calculations. By looking at the data in 

reference [12], we estimate -^m ~ 0.8 as the plateau value for large velocities, 
where the above ratio seems to converge to. 

From our model, and using the data above, we obtain 

rpm 
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which is not too far from the experimental result above. But we need to keep 
in mind that the present model is basically qualitative. It is not sophisticated 
enough in order to obtain precise quantitative agreement with the experiments. 



5 Conclusions 



A simple Langevin-like collisional model, valid for any mass ratio and inelas- 
ticities values, which assumes a bath of grains interacting with test particles, 
is used to study the steady-state average kinetic energy of a mixture of dis- 
tinct species of grains. The collisions are modeled by a Langevin-like impulsive 
term. The collisional history of the test particle is obtained and the limiting 
mean square velocity calculated. We compare the mean kinetic energy for 
each species. There results a difference in coexisting granular temperatures. 
We verify that this is due to inertial effects during a collision of massive and 
light particles. 
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